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We use spin-torque resonance to probe simultaneously and separately the dynamics of a magnetic 
domain wall and of magnetic domains in a nanostripe magnetic tunnel junction. Thanks to the 
large associated resistance variations we are able to analyze quantitatively the resonant properties 
of these single nanoscale magnetic objects. In particular, we find that the magnetic damping of both 
domains and domain walls is doubled compared to the damping value of their host magnetic layer. 

We estimate the contributions to damping arising from dipolar couplings between the different layers 
in the junction and from the intralayer spin pumping effect. We find that they cannot explain the 
large damping enhancement that we observe. We conclude that the measured increased damping is 
intrinsic to large amplitudes excitations of spatially localized modes or solitons such as vibrating or 
propagating domain walls. 


The spin torque diode effect provides an efficient 
tool to access the resonant properties of sub-micrometer 
magneto-resistive structures pQ. In particular, contrar- 
ily to conventional ferromagnetic resonance techniques, 
its sensitivity allows probing the dynamics of individ¬ 
ual magnetic solitons, such as vortices or domain walls, 
which are foreseen as information vectors in next gener¬ 
ation magnetic devices HE], By injecting a microwave 
current through the stack at resonance, spin torque can 
induce large magnetization precession of the free layer. 
Magneto-resistance then converts this precession into re¬ 
sistance variations. These resistance oscillations, multi¬ 
plied by the microwave current oscillating at the same 
frequency, give rise to a rectified dc voltage. Most spin 
torque diode studies have focused on the dynamical prop¬ 
erties of uniform magnetization configurations mm 
and few on the vibration modes of magnetic solitons cr 
:8 . Moreover, the latter studies were performed with 
metallic samples, where the magneto-resistance ratios 
are typically restricted to a few %. These low magneto¬ 
resistance ratios strongly limit the amplitude of the out¬ 
put dc signal. In the case of domain wall (DW) vibra¬ 
tions, generally confined and of limited amplitude, it is 
therefore very difficult to extract other quantitative pa¬ 
rameters than the resonance frequency of the wall [ME]. 

In this Letter, we take advantage of the large 
tunnel magneto-resistance ratios provided by MgO- 
based magnetic tunnel junctions to probe the dy¬ 
namics of a DW and its neighboring domains by 
spin torque diode effect. The junctions stack, illus¬ 
trated in Fig. [lja), Si02//buffer/PtMn (15)/Co7oFe3o 
(2.5)/Ru (0.9)/CoFeB (3)/MgO (l.l)/Ni 80 Fe 20 (5)/ Ru 


(10) (thicknesses in nanometers), gives rise to tunnel 
magneto-resistance ratios of about 30%. A synthetic an- 
tiferromagnet (SAF) prevents the formation of a domain 
wall in the CoFeB reference layer. The injection of a DW 
in the free layer is facilitated by the arc shape geometry 
of our junctions, shown in Fig. [ljb). After the saturation 
of the free layer magnetization in the y direction with a 
strong external field (300 Oe), the DW is nucleated when 
the field is decreased below 75 Oe. In Fig. [ljc) , micro- 
magnetic simulations of the full stack realized using the 
OOMMF software [151 show that the DW is first nucle¬ 
ated in the center of the wire, then displaced along the 
wire when the external field is decreased due to the stray 
field from the SAF. This scenario is in very good agree¬ 
ment with the experimental measurement of the junction 
resistance as a function of the decreasing field applied 
along the y direction shown in Fig. [ljd). The bottom 
and top black dashed lines, given for reference, corre¬ 
spond respectively to the resistance of the parallel (P) 
and anti-parallel (AP) states. Each change of the resis¬ 
tance value corresponds to a displacement of a DW in the 
Permalloy free layer. The first resistance plateau of 130 
Ft for the field range between 75 and 26 Oe (blue square 
dot in Fig. [ljd)) corresponds to the magnetic configura¬ 
tion labeled 1 in Fig. [ljc) . For lower fields, between 20 
and -35 Oe, the DW shifts to a strong and reproducible 
pinning site due to the shape at the edge of the wire m ■ 
The resulting resistance plateau at 122 Q in Fig. m 
corresponds to the magnetic configuration labeled 2 in 
Fig. [ljc). For negative fields larger than -40 Oe, the DW 
is expelled of the wire and the resistance value reaches 
the P state at 116 Cl (orange square dot labeled 3 in Fig. 
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FIG. 1. (a) Schematic of the spin torque diode measure¬ 

ment set-up including the MgO-based magnetic tunnel junc¬ 
tion with the DW in the Permalloy free layer, (b) Scanning 
electron microscope image of the sample, (c) Full-stack mi- 
cromagnetic simulations representing the different steps in the 
DW displacement, (d) Resistance versus magnetic field mea¬ 
surements. Black dashed lines: resistance values of the P and 
AP states. Black line: evolution of the resistance when a DW 
is nucleated (e.g. blue square at 26 Oe) in the middle of the 
wire (1) then pinned in the right edge (2) until its expulsion 
(e.g. orange square at -85 Oe) (3). 


0d), and magnetic configuration labeled 3 in Fig. s°». 

We probe the resonant properties of our system by spin 
torque diode in the configurations with and without a do¬ 
main wall (respectively blue square labeled 1 and orange 
square labeled 3 in Fig0 (b)). We inject a microwave 
current between the top and bottom electrodes, and 
sweep the frequency. We repeat this measurement for 
different values of applied field H y . In the absence of dc 
current, the rectified voltage V^ c has two main compo¬ 
nents (T7;: 
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FIG. 2. Normalized rectified voltage Vdc/^hf as a function 
of frequency (a) Black line: measurement at 26 Oe with the 
DW positioned in the middle of the wire corresponding to a 
resistance of 130 D. Blue line: fit with Eq.®. (b) Black 
line: measured at -85 Oe in the parallel state (DW expelled). 
Orange line: fit with Eq. §. 


1 r) 2 V r) 2 V 

Vdc =2~W + me VvWW (1) 

The first term is a purely electrical background sig¬ 
nal due to the bias dependence of the tunnel junction 
resistance, from which we extract the exact value of the 
injected microwave current in the junction. The second 
term is the spin torque diode contribution to the rectified 
voltage and arises from magnetization oscillations. The 
parameter 6 depends on the excited mode. It is related 
to the precession angle for quasi-uniform magnetization 
oscillations of the domains, and to the domain wall posi¬ 
tion for DW vibrations in a pinning potential. 

Fig. [2|a) and (b) show the rectified voltage response 
normalized by the square of the microwave current am¬ 
plitude; Vdc/I^p at two different magnetic field values. 
Fig. [2|a) corresponds to a measurement at 26 Oe, per¬ 
formed at 130 D (state T’ in Fig. [ljd)), where the DW 
is positioned in the middle of the wire. Fig. [2jb) shows a 
measurement at -85 Oe, where the magnetization of the 
free layer has a quasi-uniform configuration correspond¬ 
ing to the parallel state (state ’3’ in Fig. Qd)). The 26 
Oe response curve shows two resonance signals at 1.24 
GHz (mode 1) and 2.15 GHz (mode 2) while the -85 Oe 
response shows only one at 2.3 GHz (mode 2). We sys¬ 


tematically observed two resonance signals in the field 
range between 50 and -35 Oe where the DW is present 
whereas only the higher frequency mode (mode 2) was 
observed between -40 and -150 Oe where the DW is ex¬ 
pelled. The low frequency signal (mode 1) can therefore 
be ascribed to DW oscillations while the higher frequency 
signal can be attributed to magnetization precessions in 
the domains. When the DW is not present mode 2 is not 
strongly modified, only its resonance frequency is shifted 
by the external field applied to allow the DW expulsion. 
This means that the strongest vibrations of this mode 
are spatially located far enough from the DW not to be 
impacted by the DW vibrations when the DW is still 
present. 

Around the resonance frequency /o, the normalized re¬ 
sponse curve takes the shape of a linear combination of 
Lorentzian and anti-Lorentzian profiles [1, 14!. We fit the 
resonant signals of the mode 1 (DW) and 2 (domains) by 
the following expression: 

V*c _ A(/p - / 2 ) + -B / 2 

Ilf (/o 2 - f 2 ) 2 + (A /) 2 + j 

where / is the frequency of the microwave current 1^/ 
and the free parameters of the fit are the amplitudes A 
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and B of the anti-Lorentzian and Lorentzian profiles, a 
constant shift C, the resonance frequency /o, and the 
linewidth A. The average value of A extracted from the 
measurements by fitting with Eq. © (Fig. § at different 
magnetic fields is about 0.4 db 0.02 GHz both for mode 1 
(DW) and 2 (domains). The magnetic damping parame¬ 
ter a is related to the linewidth mm- Here, in order to 
establish the correlation between A and the correspond¬ 
ing damping a for the DW and domains resonant signals 
we perform micromagnetic simulations of the 5 nm Py 
free layer at zero external field using the OOMMF soft¬ 
ware maun. We find a resonant response for the DW 
and domains respectively at 1.7 (Fig. [3^a)) and 2.75 GHz 
(Fig. |3jb)) in good agreement with the experimental re¬ 
sults. To reconstruct the spin diode signal we sweep the 
frequency of the microwave current Ihf between 0.9 and 3 
GHz and extract the resulting magnetization oscillations. 
We plot as a function of the frequency the average am¬ 
plitude of the product between the x component of the 
magnetization oscillations M x and lhf= Io*sin( 27 r/£ + (j>) 
(Fig. |3|. The </> parameter (respectively equal to ^ and 
0 in Figs. j3^a) and (b)) is related to the resonant signal 
shape depending on the symmetry of the exciting force 
(Slonczewski or Field like torque [l]) and does not influ¬ 
ence the linewidth. For a damping parameter <a=0.01, 
the linewidth extracted from the simulated spin diode 
signal by fitting with Eq(2] is 0.15 ± 0.002 GHz both for 
mode 1 and 2 (Fig. [3|. 

For the particular case of the DW resonance, in a ID 
model, the linewidth A is related to the damping param¬ 
eter a by: 


A 7o aH d 

27r 


( 3 ) 


> 

35 0 

^ 40 

CD 

"co 30 

CD 

O 20 

10 

Id-10 

03 


yj 1.0 1.5 2.0 2.4 2.6 2.8 3.0 

Frequency (GHz) 



FIG. 3. Spin diode signal reconstructed from micromagnetic 
simulations of the Py free layer at zero external field, (a) 
Black line: Reconstructed spin diode signal in the 0.9-2.2 
GHz frequency range, showing the resonant response of the 
DW (mode 1). Red line: fit with Eq. ([2]). (b) Black line: 
Reconstructed spin diode signal in the 2.4-3 GHz frequency 
range showing the domains resonant response (mode 2). Red 
line: fit with Eq. ©• Insets: Spatial distribution of the am¬ 
plitude of the magnetization oscillations dynamic component 
from the SpinFlow3D mode solver associated with the two 
observed modes. The color scale blue-green-yellow-red corre¬ 
sponds to the amplitude of the mode. The white arrows show 
the direction of the magnetization. 


to the resonant response of the mode 2. For the 
mode 1, the resonant response is associated to a transla¬ 
tional vibration mode of the DW (inset in Fig. [3|a)) (19: . 


where 70 is the gyromagnetic ratio, A the linewidth 
in Hz, and H d the demagnetizing field, that we calculate 
to be equal to 0.45 MA/m in our samples where the 
magnetization of the free NiFe layer is 0.47 MA/m. 
Based on this ID model (Eq. [ 3 ]), the expected linewidth 
from a damping of 0.01 is 0.15, in good agreement with 
the linewidth extracted from micromagnetic simulations. 
This result indicates that the ID model is relevant to 
evaluate the damping extracted from the experimental 
linewidth measurements for the DW. Moreover, the 
fact that the ID model also allows to describe the 
domains damping hints once again to the fact that the 
magnetization oscillations in the domains are spatially 
localized and confined. This is confirmed by the spatial 
distribution of the modes calculated using the mode 
solver from the SpinFlow3D package na uni where 
the full stack is considered [21], shown in the insets 
of Fig. §a) and (b). Indeed, the most important 
dynamic component of magnetization in the domains 
is confined and located in the left edge of the wire 
(inset in Fig. |3jb)). This edge mode is associated 


Based on the previous micromagnetic simulations and 
the ID model expression (Eq. [ 3 ]), the linewidth extracted 
from the measurements corresponds to a damping a of 
0.026 ± 0.001 both for DW and domains magnetization 
oscillations. This value of a is more than twice the value 
that we have measured by FMR on the unpatterned 
film (0.01T 0.0015) and that is typically observed in Py 
films (0.01 [TT]) (see supplementary material). For each 
of the five samples we have measured we have obtained 
similar large damping values for the vibration modes, in 
the range 0.019 to 0.028. 

Several mechanisms can be involved in this large in¬ 
crease of the effective damping compared to classical 
FMR measurements on thin films. We consider possi¬ 
ble effects at stake and estimate their relative contri¬ 
butions. Surface roughnesses or edge defects are often 
pointed out [22] but it has been recently shown by mi¬ 
cromagnetic simulations that the resulting damping is an 
order of magnitude smaller than what we observe [23, [24] . 

Another potential source of strong damping increase, 
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especially in the case of a DW, has been pointed out by 
Ndjaka et al. [25]. Indeed, in trilayer structures, the stray 
field generated by DWs in one ferromagnetic layer dipo- 
larly couples to the other ferromagnetic layer, leading, 
in the case of very strong couplings, to DW and domain 
duplication [23 EZ]. In the case of weak coupling only 
a local perturbation of the magnetization in the refer¬ 
ence layer is observed [25]. This is what we expect in our 
magnetic tunnel junctions as the SAF prevents domains 
duplication. This magnetic shadow in the reference layer 
coupled with the DW dynamic in the upper one results in 
an increased damping. To confirm and quantify this be¬ 
havior, we have performed micromagnetic simulations of 
the full stack [21] using the OOMMF software [15]. Figs. 
[4ja) and (b) show respectively the magnetization of the 
Py free and the CoFeB reference layer. As expected, we 
observe that the magnetization of the reference layer re¬ 
veals a non-uniformity located just below the DW of the 
free layer. We extract the effective DW damping anw 
from the DW damped oscillations (inset of Fig. S c)) as 
a function of <a re /, the damping of the SAF layers (Fig. 
[4]^c)) (see supplementary material). 

In the simulations, we set the damping of the free layer 
to 0.01, and vary the damping of the SAF layers between 
0 and 0.1. Figure [4] (c) shows the dependence of oldw 
on a re f for two values of the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interlayer exchange coupling in the SAF: 
—1.10 -4 (filled symbols) and — 3.10 -4 (open symbols) 
J/m 2 and for three values of the reference CoFeB layer 
magnetization: 0.82 MA/m, 1.03 MA/m and 1.2 MA/m. 
In each case, we observe an increase of with <a re j, 

thus confirming the dynamic coupling between the DW 
and its magnetic shadow. However, in our experiments, 
the typical values of the RKKY exchange coupling, a re f 
and Ms for the CoFeB reference layer are respectively 
-0.1 erg/cm 2 , 0.01 and 1.2 MA/m [28] [29] . As shown 
in Figure[4] (c), the corresponding oldw is about 0.0115, 
meaning that the DW damping increase due to dipolar 
coupling is one order of magnitude smaller than what we 
observe. 

An alternative phenomena that can lead to strong 
damping increases is intralayer transverse spin pumping 
[23] m [30] EB. In a conducting ferromagnet, the 
conduction electrons carry away the excess angular 
momentum of a precessing magnetization. For spatially 
inhomogeneous magnetization dynamics, a nonuniform 
spin current is induced, resulting in a spatial dependance 
of the dissipative flow of conduction electrons in the 
magnetic layer itself. This can give rise to an enhanced 
damping in isolated magnetic nanostructures. This 
intralayer spin current results in the following additional 
nonlocal torque in the Landau-Lifshitz-Gilbert equation 
[24 : gt(j& x where gt is the transverse spin 

conductivity. In Ref. [30], the increase of damping 
due to intralayer spin pumping is estimated by ^(^-) 2 , 



FIG. 4. OOMMF micromagnetic simulations of the Permalloy 
free layer with the DW (a) and the CoFeB reference layer (b). 
(c) Dependence of the DW effective damping on the reference 
layers damping. The filled (resp. open) symbols correspond 
to an RKKY exchange coupling in the SAF of —1.10 -4 J/m 2 
(resp. —3.10 _4 J/m 2 ). The value of the CoFeB reference layer 
magnetization is varied between 0.82 MA/m, 1.03 MA/m and 
1.2 MA/m. Inset: The black solid line is the M x component 
of the NiFe layer, describing the damped oscillations of the 
DW. The red solid line is the fit (see supplementary material). 


where 77 = , is a material dependent parameter 

related to gt of Ref. [24] . The parameter A is the 
wavelength for the magnetization pattern we consider. 
This means that the more magnetization dynamics is 
spatially located the smaller is A, and the more the 
increase of damping is important. The spatial profile of 
the calculated modes shown in inset of Fig. [3] can give 
us a rough estimate for A, typically 100 nm for both 
DW and domains. By taking the parameter 77 equal to 
0.76 nm 2 for Permalloy with Ms=0A7 MA/m and the 
conductivity Go = ( bfiVtcm ) _1 [30], we find that the 
associated damping increase is 0.003. Recent analytical 
calculations confirm that for transverse domain walls 
damping increases due to intralayer spin pumping are 
very small [32] . 

Finally, spin transfer torque can give rise to large am¬ 
plitude motion of magnetic objects. In our experiments, 
the injected microwave current density is 1.43 10 5 A/cm 2 , 
about 10 times smaller than the critical current densi¬ 
ties needed to depin the domain wall m- From the 
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measured rectified voltages in Fig. [2] we find that the 
amplitude of domain wall translation at resonance is 4 
nm, a value far from negligible compared to the domain 
wall width of 100 nm. It has been shown experimentally 
that non-linear contributions, that are not accounted for 
in micromagnetic simulations, can give rise to very large 
damping increases in nanostructures [321 • By elimination 
it seems that only these non-linear contributions appear¬ 
ing for wide vibration amplitudes can be responsible for 
the very large damping values that we measure. Our 
results are in phase with the high damping values gen¬ 
erally derived from magnetic field or spin torque driven 
domain wall motion over large distances [34]. They also 
show that the damping of simple objects like transverse 
domain wall in a standard material such as Permalloy are 
still far from being understood. 

As a conclusion, we have shown that by combining the 
spin diode effect with the large magneto-resistance ratios 
provided by magnetic tunnel junctions, we can quan¬ 
titatively analyze the resonance of individual magnetic 
confined modes. The effective damping of the domain 
wall and edge modes extracted from the rectified spin 
diode signal is more than twice increased compared to 
the one extracted from FMR measurements on the ex¬ 
tended stack. We ascribe our observations to non-linear 
contributions to the damping appearing for large ampli¬ 
tudes of vibration. Our results underline that damping 
mechanisms are still far from being elucidated in simple 
systems. They also show that non-linear contributions 
should be taken into account when modeling large am¬ 
plitude vibrations or translations of solitons in magnetic 
nanostructures, in particular for memory and logic appli¬ 
cations. 
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